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Abstract. - We obtain the memory kernel of the generalized Langevin equation, describing a 
particle interacting with longitudinal phonons in a liquid. The kernel is obtained analytically 
at T — Q°K and numerically at T > 0°K. We find that it shows some non-trivial structural 
features like negative correlations for some range of time separations. The system is shown to 
have three characteristic time scales, that control the shape of the kernel, and the transition 
between quadratic and linear behavior of the mean squared distance (MSD). Although the 
derivation of the structure in the memory kernel is obtained within a specific dynamical model, 
the phenomenon is shown to be quite generic. 



The Generalized Langevin equation is a powerful tool for the study of dynamic properties 
of many interesting physical systems. It is widely believed that the memory kernel is some 
decaying function of time, with no interesting features, except for the characteristic decay 
time. An example is the simple Gaussian form used in the literature [Q, ^]. Some molecular 
dynamics simulations yield, on the other hand, non trivial structure of the memory kernel ||. 
A result of the same nature was obtained analytically by Chow and Hermans p| . Their memory 
kernel is positive only at t=0 and negative elsewhere. The derivation assumes, however, that 
the Brownian particle interacts with a viscous fluid and therefore it is not a first principle 
derivation. Furthermore the final form of the kernel obviously violates the requirement that 
its Fourier transform be positive for all ui. 

The history of attempts to actually derive the memory kernel of a generalized Langevin 
equation from a first principle underlying Lagrangian or Hamiltonian description starts with 
the paper of Feynman and Vernon |^|. Feynman and Vernon couple the Brownian particle 
linearly to a system of harmonic oscillators, obtaining an exactly solvable problem. Then the 
memory kernel is derived from the microscopic parameters. As explained by Schwartz and 
Brustein [Q, it is difficult to envisage a situation where the Brownian particle is not bound to 
some small region in space, where a linear coupling could be justified. A number of authors 
followed in this direction. Zwanzig Q and Lindenberg || developed general formalisms to 
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obtain the memory kernel, but their actual implementation, in obtaining the memory kernel 
in terms of the parameters of the Hamiltonian, is restricted to systems with linear coupling 
of the Brownian particle to the degrees of freedom of the thermal bath. The formalism of 
Zwanzig, combined with the work of Mori [ll| , has been extensively used in the research 
that followed (for examples see: Goodyear and Stratt p2[ , Guenza jlS, Heppe |I3 and lee 
III)- 

It is clear that a microscopic derivation of the memory kernel is needed, for the cases where 
the particle is not bound, to determine whether and under what conditions the kernel has some 
interesting structure. In this paper we derive the memory kernel for a particle interacting in a 
realistic way with the longitudinal density waves of the system in which it is immersed. (The 
interaction has the same form as the electron-phonon interaction). We find that the structure 
of the memory kernel has interesting features. The most interesting is the existence of a region 
of time for which the memory kernel is negative. Although our result is obtained for a specific 
model, we show later that it must be generic. 

Former work by Munakata considers a Lagrangian describing an impurity interacting with 
an elastic periodic lattice and obtains from it the shape of the memory kernel fl6[| . His 
derivation is based, however, on the long time form of the MSD, that is linear in time. To 
refine his result to be also valid at very short times, one needs to obtain the full time dependence 
of the MSD. In this article we present a method that yields the mean squared distance and the 
memory kernel, both over the full time range. Our starting point is the Lagrangian describing 
the interaction of a particle with longitudinal phonons in a liquid; considered in refs. |7], 17 



We consider a system of a particle immersed in an infinite idealized liquid and interacting 
with it's longitudinal phonons through a general two body interaction, u(q). 
The Lagrangian of the system is [jl7| : 

■ ~ lv(q)p(q)p(-q) KO) [P (-0) ■ /(-$)] ) (1) 

1 .2 

dq(u{q)p{—q)exp(iq ■ x)) + -M x 

where m is the mass of the liquid particles, ~p is its average number density, J(q) and p{q) are 
the Fourier transforms of the current density and the number density respectively, the Fourier 
transform (FT) of the two body (effective) potential between particles of the liquid is V(q) and 
the FT of the interaction between the Brownian particle and the liquid is u(q). The Lagrange 
multiplier fj,(q) is introduced to impose the eq. of continuity, x is the coordinate of the particle 
and M is its mass. 

Assuming that in the absence of the particle the distribution of the degrees of freedom of 
the liquid is given by a Gibbs distribution at temperature T, Brustein et al. found Jj"7| that 
the particle obeys a generalized Langevin eq. . 

M S(t) = - [ dt'~f(t - t') x (<') + F(t) (2) 

J — oo 

The average of the random force F is zero, and the force - force correlation is related to the 
memory kernel. 



F) = (3) 
(F(t)-F(t'))=3K B T 7 (t-t') (4) 
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The memory kernel was also given in terms of the interactions, phonon frequencies, u>(q) and 
the MSD. 

lit -0 = i/ d ^V^ C0S [«(«)(* - 01 e-*« A '(*-0) a ) (5) 

In ref. |l7j a numerical method for calculating the kernel was suggested. The numerical 
method was based on iterations. To calculate the MSD one had to evaluate four dimensional 
integrals, involving Laplace transforms, a difficult numerical procedure. In this article, we 
suggest a slightly different approach that involves only one dimensional integrations 

The first step is to obtain the MSD from a generalized Langevin eq. with a given memory 
kernel 7. To solve eq.(2), we transform it to Fourier space. We do it by inserting a step 
function: 

/>oo 

Mx + ~{{t-t')Q{t-t')'x{t')dt' = F(t) (6) 



Solving directly for x (u>), we obtain: 

(M*) a )= -If^fcfc/ <>H->K))(^D(e^-l) > (?) 
X ' 2ttJ_ oc J_ oc (y/2^iujT(uj) - Muj 2 )(V2^iuj>r(uj') - Mlo' 2 ) 

A 

where T(oj) is the Fourier transform of j(t)Q(— t), and F (00) is the Fourier transform of F(t). 
Using the fluctuation dissipation relation (4): 

A A A 



(F (lo)- F (—u)')) = 3KbT\/2tt 7 (lo)S(lu — lj'), a straightforward calculation yields 
2X 6K B T f°° ^ 7(w)(l-co«(w*)) 



(Ax(t) 2 ) = — / dw—~ — ^^ / ; ' -, (8) 

A 

where 7 (uj) is the Fourier transform of r f(t). 

Rescaling the integral on the right hand side of the above by defining y = u>t, it is clear that 
the long time dependence of (Ax(t) 2 } is linear in t. The short time dependence is quadratic 
in t as expected and given by 



3K B T f°° , 7 (w) 



2ir J^oo (Mlu) 2 + 2ttT(uj)T{-lj) 



ipt 2 (9) 



The two required quantities 7(i) and (Ax(<) 2 ) are obtained now by an iterative procedure 
from the two coupled eqs. (5) and (8). In order to be specific we need, however, to determine 
the potentials u(q) and V(q). We take for V(q) a constant, V(0), corresponding to a S function 
in real space. For u(q), the potential between the immersed particle and the fluid particles we 
take a potential that has a finite range of the order of the size of the particle. We choose a 
Gaussian form u(0)e~ a q . The choice of V(q) yields for u)(q), the phonon spectrum 



u(q) = \ ^P-\q\=c\q\ (10) 
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Fig. 1. - The memory kernel in dimensionless units ( 70 = v*(q=o)a 5 •) as a function of 
dimensionless time, for fixed values of T\ and T2 . 



where c is the sound velocity. 
The iteration procedure is defined by the following eqs. 



G(")(q,^) = e -^ 2 « A ^» 2 > ( ' ,, (11) 

7 (n) (t) = A _lJ dqq 4«M cos[ujm)]G (n) {q ^ (12) 

1 r° 

r( n \uj) = -= ^ n \t)e- i0Jt dt (13) 

V 27T J-00 

7™ } (w) = / 7 ( ™ ) (t) e -^ t dt (14) 

V27T J 

and finally 



(A4<)) - loo ^ ((M c)2 + 2^r(«)Mr(«)(-c)) (15) 



(n) 



For the specific potentials we chose, the memory kernel is: 

lit) = \y^j coa^e-W^ (16) 

Calculating the integral while paying attention to the fact that (Ax(t) 2 ) does not depend on 
q ,we obtain: 



nu 2 (0) ^ {\2A 2 - \2B 2 A + B 4 ) e 



where: 



A = 2a 2 + 1 (Ax(t) 2 ) 
B = ct 

Since the MSD vanishes with temperature, the calculation yields an explicit analytic expression 
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Fig. 2. - Dependence of the memory kernel on ti and time for small t\. 
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Fig. 3. - Dependence of the memory kernel on r 2 and time for large t\ . The fig. presents three 
lines that almost coincide. 



for j(t) at T = 0°K, with A — 2a 2 . In ref. jlTj it was found that the friction coefficient, 7, 
vanishes at T — Q°K. It is interesting to note that indeed 



j(t,T= 0)dt = 



(18) 



The results at finite temperatures can be described in terms of three natural time scales: to, 
t\, and ti- The time to is the time needed by sound to traverse the distance a, that is the 
size of the particle; to = -. The time t\ is also related to the size of the particle and is 
proportional to the mean time that takes the particle to move its own size in the short time 
behavior regime. It is given by t\ = (eq. 8). The third characteristic time ti is given by 

t2 = — , where 7 =7 (u> — 0) (Notice that the friction coefficient of the regular Langevin eq. is 

■^E times 7). In fig.(l) we depict a typical kernel that exhibits non-trivial structure. A fast 
decay is followed by a negative region, a small positive peak and then a final decay. The most 
interesting feature, in the structure of the memory kernel, is the negative region, of course. 
The shape of the memory kernel depends on the two dimensionless parameters t\ = |j and 
t 2 = t 2 -. Fig. (2) gives the dependence of the kernel on ti for fixed t\ « 1. We see a very 
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weak dependence of the structure on t 2 ranging from much below to much above unity, without 
any particular interesting features. When t\ is kept fixed at values above 1 and r 2 is varied, 
we see (Fig. (3)) that although the structure of the memory kernel is very different from the 
structure for n << 1 the dependence on t 2 is still weak. In fig. (4) on the other hand, we see 
the dependence of the kernel on t\ with r 2 fixed. The striking feature is the transition from 
a simple decay at small t\ to an interesting structure with a negative region that becomes 
more pronounced as t\ is increased. The origin of the negative region can be traced back 
to the discussion in refs. |?J and |l7) that shows that at zero temperature / °° j(t,T = 0)dt 
vanishes. Therefore, it is clear that at zero temperature j(t) must have a negative region to 
compensate for the contribution of the positive region. It is clear therefore, that for low enough 
temperature we must also have a negative region. We are dealing in this paper with a specific 
model describing a particle interacting with longitudinal phonons in an idealized liquid. We 
would like to stress at this point that the fact that the memory kernel has a negative region 
at low temperature is generic and not just an artifact of the specific model we consider. The 
reason is that the fact that at zero temperature J Q j(t,T = Q)dt = is independent of 
the model and quite general. Therefore, the above conclusions concerning the existence of 
a negative region apply to the general case. In fact, this result is supported by molecular 
dynamics simulations that show similar behavior of the memory kernel ^| . Returning now 
to our model, we have to define low enough temperature in term of a dimensionless quantity. 

A careful inspection reveals that the relevant dimensionless parameter is and it is easy 

_ i 

to show that t\ — d ■ ( ^/JT ) 2 ■ Indeed d is a function of n and t 2 , but it does not vanish 
or diverge as t± or t 2 tend to infinity. It is obvious now why when t\ is increased a negative 
region develops in the kernel. The parameter t 2 , on the other hand, can be increased without 
increasing t\ (or equivalently, decreasing \fji ) ■ This explains the weak dependence on t 2 
when T\ is fixed. The physical origin of the negative region is the finite size of the particle and 
as result, the interaction of the particle at one point with phonons emitted earlier at another 
point. (The emission of phonons is the mechanism by which the particle looses energy. The 
interaction with phonons can contribute to the energy of the particle). It is interesting to work 




Fig. 4. - Dependence of the memory kernel on t\ and time. 



out the times t where the kernel vanishes. This will tell us also when a negative region exists. 
We will demand that "f(t) given by eq.(17) vanishes and will assume that ^(t) vanishes when t 
is small enough, so that the expression to be used for the MSD is Ax 2 (t) — tpt 2 (eq.9). Using 



G. FRENKEL et al.\ THE STRUCTURE OF LANGEVIN'S KERNEL ETC. 



7 



x 10 4 




-40 -30 -20 -10 10 20 30 40 



Fig. 5. - Typical MSD, presented in dimensionless units (a is the range of interaction between 
a particle of the medium and the Brownian particle). 

the short time form of the MSD is reasonable since at zero temperature the MSD vanishes, and 
therefore we are interested at finite temperature in the region where the MSD is still small. 
Solving the system of eqs. 

12A 2 - 12B 2 A + B 4 = 

A w 2a 2 + \<pt 2 
6 

B = ct 

we find that the kernel vanishes at two points, one point, or none at all; depending on t\. 

*(7 = 0) = J 2 2 t, 7- = J 2 2 1 j (19) 

V 12±V96 6^ V 12±V96 1 6 

When n >> 1 a negative region will exist; however, when n << 1, the negative region will 
disappear. In fig. (4), curves (l)-(3) correspond to large n and indeed the memory kernel 
vanishes at two points; the calculation shows that curves (4) and (5) should have only one 
zero point as can be seen in the fig.; and curve (6) should have no negative region at all. 

In fig. (5) we present a typical dependence of the MSD on time. As expected, the MSD is 
quadratic for short times, and linear for long times. It is easy to show, from eq.(8), that the 
linear long time dependence is (Ax(t) 2 ) = — 12J f sT — • t, and we have shown the short time 

v/2¥-7(w=0) 

dependence to be (Ax(t) 2 ) = <pt 2 ; therefore we can easily show that the transition between 
them occurs at a typical time scale which depends only on the three known time scales. When 
T\ < 1 the transition time is roughly t 2 . 

Our calculation is classical and clearly at zero temperature quantum mechanics prevails. 
Nevertheless it is obvious that can be made small enough while the system still is not 

quantum mechanical. Therefore, this generic effect should be experimentally observable. 
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